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ABSTRACT 


This document summarizes the mathematical model and test results 
developed for the Vertically Reciprocating Auger (VRA). The VRA is a 
device capable of transporting cuttings that result from below surface 
drilling. It was developed chiefly for the lunar surface, where 
conventional fluid flushing while drilling would not be practical. 

The VRA used only reciprocating motion and transports material 
through reflections with the surface above. Particles are reflected 
forward and land ahead of radially placed fences, which prevent the 
particles from rolling back down the auger. Three input wave forms are 
considered to drive the auger. A modified sawtooth wave form was 
chosen for testing, over a modified square wave or sine wave, due to its 
simplicity and effectiveness. 

The three dimensional mathematical model predicted a sand 
throughput rate of 0.2667 pounds/stroke, while the actual test setup 
transported 0.075 pounds/stroke. Based on this result, a correction factor 
of 0.281 is suggested for a modified sawtooth input. 
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PROBLEM STATEMENT 

The objective of this project is to develop a mathematical model for 
a mechanism capable of transporting drill cuttings up a vertical shaft 
using only linear reciprocating motion. 

Background 

It is anticipated that there will be a need to drill into the lunar 
surface. It is also anticipated that conventional methods of cutting 
removal (fluid flushing) will not be practical. Therefore a method of 
cutting removal that relies on mechanical dynamics is needed. To 
accomplish this task, an accurate math model and prototype must be 
developed and tested. 

Performance Objectives 

The objective of the math model is to accurately predict the 
amount of material the auger is capable of transporting. The math model 
must be constructed to allow for variations of parameters such as pitch, 
fence height, fence separation, gravity and auger velocities. 

The prototype must be capable of transporting material and 
providing information that can be used to improve the math model. 

Constraints 

The mechanism is allowed to move linearly along its center-line 


axis. 



MATHEMATICAL MODEL 


Introduction 

Th© vertical reciprocating auger (VRA) is able to transport material 
by virtue of its geometry and motion. The auger studied is a common 
helical screw type with one pitch per turn and pitch equal to outside 
diameter . However, the relations that follow are generalized to predict 
particle dynamics of any screw type auger, through a choice of the 
appropriate geometric parameters. In this study, the auger reciprocates 
in the vertical axial direction only, coincident with a gravitational force, 
and is not subjected to a torque about its axis. The VRA differs from 
normal screw type augers by having an outside fixed sleeve and radially 
located vertical fences attached to the upper side of the helix. The fences 
prevent particles from rolling down the auger during transport. 

To derive the mathematical model of particle dynamics, the motion 
of a single particle is studied throughout a single cycle of input motion to 
the auger. Three different reciprocating inputs are proposed; namely 
sinusoidal, a modified square wave and a modified sawtooth. Next, 
through multiple cycles of a reciprocating input, a transport rate of this 
particle up the auger will be predicted. The analysis will then be 
extended to include many particles, with a correction factor suggested 
from actual empirical results. This single particle analysis is reasonable, 
as the math model will show that particles at different radial locations 
along the auger will trace different paths in 3-space. 

The particle motion is first analyzed by "unwrapping" a turn of the 
auger for a linear, or 2-dimensional representation (Fig 1) . The "floor" 
and "ceiling" of a single auger turn make an angle, with the horizontal, of: 

O = tan ' 1 ( P / EnRj ) 0 ) 

where : P = Auger Pitch 

Ri = Any radius, i , from the central axis. 



A particle is able to advance to higher fences through reflection off 
the ceiling. A particle begins at rest ahead of the first radial fence, and is 
accelerated upward with the auger. As the auger reaches a steady 
speed and begins to decelerate, the particle separates from the surface 
and moves upward. With sufficient kinetic energy, the particle impacts 
the ceiling at an incident angle 6j to the ceiling normal and is reflected at 
an equal angle, 0 r , to the normal. The total angle of reflection from the 
incident path is p = 20. The particle now has a horizontal component of 
Vsin p, which will allow it to clear the next fence up the auger. This 
process is repeated and material is transported vertically up the auger. 

Inputs 

A. Sinusoidal 

The first reciprocating input studied was the sine function (Fig 2), 
in which the particle will reflect once per cycle. Figure (3) shows the 
particle position at tn time increments as follows: 

to: Auger and particle are at bottom of motion, velocity = 0. 

ti : Auger and particle now at half of amplitude. Auger begins to 
decelerate. Particle separates, travels upward at velocity vi p , 
which is the maximum velocity of the auger, occurring at y=0. 

t 2 : Auger reaches top of amplitude, begins to accelerate 
downward, while particle continues upward at approximately 
vip. 

t 3 : At the ideal critical frequency of this system, the particle 
traveling upward at vip impacts the auger, traveling at (-vi p ) 
at time t 3 . The particle now reflects off of the ceiling at 2vi p . 

t 4 : As auger is decelerating in its downward motion particle 
impacts bottom, loses energy and settles above the nest 
fence. 

t5: Auger and particle are at the bottom of travel - same as to. 



The sine input is advantageous since a simplistic driving 
mechanism could be constructed from a rotary driver. However, this 
input was not chosen since "tuning" the mechanism to a critical frequency 
would very difficult. 

B. Modified Square Wave 

The next input, that was considered, modeled a double-acting 
pneumatic or hydraulic actuator. This system input takes the form of a 
modified square wave (Fig. 4 ), and involves an upward and a downward 
throw. With this input, the particle can complete two reflections per cycle 
(Fig 5 ) as the following time increments describe: 

to: Auger and particle are at bottom of motion, velocity = 0 . 

t-j: Auger is nearing top of motion, and beginning to decelerate -- 
particle separates. 

t 2 : Auger has been brought to stop -- moving particle strikes 
ceiling of auger level at vi and is reflected off at approximately 
vi. 

t 3 : Particle impacts floor of auger, losing energy, just above the 
next fence and settles. Auger still at rest at its uppermost 
position. t3 --> t4 can be any length of time. 

t4: Auger is rapidly thrust downward, leaving particle suspended. 
Throw must be slightly greater than pitch to account for small 
drop of particle. 

t 5 : Ceiling of auger strikes particle and particle is reflected ahead 
of the next fence. 

t 6 : Auger is at rest at bottom and particle impacts, loses energy, 
and settles. 

ty: Cycle is completed--same as to- 

This square wave input has several advantages over a sinusoidal 
input. First, hydraulic and pneumatic actuators are simple, reliable, and 
easy to control. Also, the motion is no longer frequency dependent; in 
fact, any length of time can occur between t3 and t4 or t 6 and t7. Finally, 
two jumps are now obtained for every cycle of actuation. 



The main disadvantage of this input is the length of stroke, which 
must be slightly greater than the pitch to impact the suspended particle. If 
the stroke is too small, so that the suspended particle does not strike the 
ceiling, it may bounce back over the fence and defeat the entire stroke. 

C. Modified Sawtooth 

The final input studied, and used for the actual test, is the modified 
sawtooth wave (Fig 6). This wave describes the output of a single acting 
pneumatic or hydraulic actuator. This sawtooth wave works similarly to 
the square wave for the upward throw, but decays with the exhaust of the 
actuator, rather than a downward thrust. The particle will now complete 
one reflection per cycle (Fig 7), as the following time increments 
describe: 

to: The instant of time, just before the auger is accelerated 
upward, when the particle and auger are still at rest. 

t-,; Auger is nearing top of motion, and beginning to decelerate— 
particle separates. 

t 2 : Auger has been brought to stop-moving particle strikes 
ceiling of auger level at vi and is reflected off at approximately 
vi. 

t3: Particle impacts floor of auger, which is moving downward, 
and loses energy. Between t3 and t4 the particle settles ahead 
of next gate. 

t4: Particle and auger come to stop at bottom of motion. 

ts: Cycle is completed-same as to- 

Although this input only moves the particle one jump per cycle, the 
sawtooth wave was chosen for testing due to its effectiveness and 
simplicity. Like the square wave, the sawtooth wave is not frequency 
dependent, and the efficiency of mass transport per cycle is not affected 
by frequency. The sawtooth wave can use a shorter stroke, since it only 
needs to transfer enough kinetic energy to particles to cause a reflection. 
At the lower limit, an infinitesimally small stroke can be used with a high 


impulse of energy to the auger. In actuality, the stroke length varies with 
the force and acceleration of the actuator. 

Two Dimensional Dynamic Analysis 

With the auger "unwrapped" at a particular radius Rj, a dynamic 
analysis can be performed on a single particle to predict its motion. A 
force balance is first done on the auger/particle system for the vertical 
throw at the actuator, from rest at the bottom to point of separation. 

m SyS Xfiy-fe = — roSy-fi? “ drag ( 2 ) 

where m sys includes the mass of the particles and auger, and Fact is the 
actuator force. Drag on the particle is assumed negligible. Rearranging, 

g* 9 v ' 

where x S ys is constant throughout the throw. 

Assuming that the particle does not separate until the actuator throw, x, is 
completed, the relationship for x is : 

X 

f % f 

X = x 0 + + .5Xf yf -k (4) 

where xo = position at bottom = 0, and vo — velocity at bottom = 0. 
Actuator throw time, t, is : 

t = y~8x f *Sy* ^ 


and separation velocity, vl , is: 



The particle now separates and is undergoing freefall in a vacuum with 
an initial velocity, vl . The distance to impact the ceiling is the pitch, P, 

Where: P . - .5,t» 


or, 


5gt^ - V, t 


P = O 


( 7 ) 


solving for the time of impact, t, with the ceiling : 

t = v u ± /~<v L >* - 2 gP (0) 

g 

and the impact velocity, v-i , with the ceiling is : 

v, = (9) 

The particle is assumed to impact perfectly elastically, and recover 
the entire vi velocity, now at an angle p = 20 in a direction defined as Ni 
(Fig 8). A direction N 2 is defined orthogonal to Ni, to form a rotated set of 
coordinate axes. The true y--axis is defined positive in a direction 
opposite to gravity, and the x--axis is defined horizontally. 

The Ni and N 2 differential equations of motion are: 

l\i * = g ( c o > ("1 0) 

Ng =-g(sin(?) (11) 

with initial conditions 

N, ( t 0 > = 0 

N , ( t o ) = V , 

N*<t 0 > = 0 (12) 

Ne < t 0 > = 0 


with to being the time immediately after impact. 



Solving thG Ni and N 2 differential equations of motion subj©ct to th© 
initial conditions : 


g ( c o sft ) 
E 

(13) 

-g i 5 i ) t * 

(14) 


£ 


Now converting Ni and N 2 positions to x--y coordinates at the at the 
ceiling impact point yields : 


X 


y 




N, ( b i n <( ) 


Nj?( cosft ) 

(15) 

( s i nft ) 

(16) 


where :y = 90 - p 

A computer program of x--y position was plotted (Appendix) for 
expected values of actuating force. The velocity is high enough, and the 
corresponding travel time is short enough, that gravity has a negligible 
effect on the particle path. The path appeared to be in the Ni direction, 
or in other words, a straight, as opposed to a parabolic, path seemed to 
be traced. 

An assumption can now be made that the reflected paths can be 
approximated by straight lines for the the 3--dimensional analysis that 
follows, when the auger is "rewrapped". 

3-Dimensional Position Analysis 

Since the actual path of a particle does not lie in a single plane, a 
3-dimensional analysis is required. Figure 9 shows a top view of an 
auger segment with a path projection of the particle, starting at radius Rj, 
at the i(th) position along a fence. 



Before the particle can impact the floor of the auger after a ceiling 
reflection, it impacts the outside sleeve. Viewed from this same top 
figure, the particle reflects about the surface normal at an angle of 2(a). 
Alpha varies from zero degrees to ninety degrees as the starting point, Rj, 
varies from zero to Ro- Orientations of planes A:A and B:B are shown on 
the top view, normal to the incident path and sleeve reflection path, 
respectively. 

View A:A is shown on figure 10 and is identical to the 2- 
dimensional model view, except that the sleeve impact point, I, is shown, 
along with a projection of reflected path (dashed line). 

From the top view, 

d, = -/~ : +- 'r; (17) 


and from A:A : 

F , = D j / si nft 


= F,COSft 


y <r„ )* + <r; >* (18) 

t an£ 


but, 


tancx = R; / V ( R,, )«? <R*. 


and, 


p 5 = r; / tancx-tan^ 


(19) 


View B:B is also shown on Figure 10 , and is taken normal to the 
path between I and floor impact. The Z direction remains the same, and 
P-i is repeated. It is noted that the same p angle between the vertical 
climb and ceiling reflection is also between the Z axis and the sleeve 
reflection path. Ci is defined as : 

c, = p - pi (20) 

= m, + 


where r r = ramp rise from lift-off 
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Define Z p = r r . Z p can be solved from the intersection of two 
relationships, as shown in figure 11 . The straight line in this figure 
represents the particle path from view B:B. The parabolic curve 
represents the climb in the Z - direction of the auger for position along 
Ybb- For the first relationship, 


Zp = C, + my^ 


(21) 


where m = slope 

Y(|) = a distance along Ybb as <|> is varied 


but since, 


O = . 5f5 = tan-i<P / PttR; > (22) 


and, 


= tan(J = tan<Stan-’CP / EttR; > > ® -1 / m (23) 

A Zbb 

so that, 

-y* 

+ P - R; / tsnn-tanft ^24) 

tan ( 2tan"l ♦ i P / E'ttR ', ) > 

The second relationship can be derived with the aid of figure 12 . 
This figure is a top view of the auger with the horizontal axis going in the 
Ybb direction, starting at the sleeve impact point, I. <D is a variable 
representing the angle between the fence reference line and any other 
radial line. By varying $ from (90 - a ), which is at the sleeve impact 
normal, to 2k, a relationship between Ybb a nd $ can be found, r is 
defined such that : 

T = f + ZO - 90° 



12 


and from the figure's geometry : 

y# = R<>coso - R| / tan! t ) (25) 

and since a = sin-i <r; / n a ) • 


then 

= R 0 cos < s i n~ 1 ( R; / R e ) ) - 


R; 


(26) 


tan ( $ + 8 < S i n~ 1 R \ / R„ ) - 90°) 

However, for a constant Rj. the rise of the auger, Z r , varies linearly with 
angle of twist, 4> : 


2 r = P$ / 360° 


(27) 


or by rearranging, 

$ = 360° 2 r / P 


(28) 


y+ 


which can be substituted into eqn. (26) to yield : 

R© c o s < s i n"‘ 1 ( R * / R 0 ) - 

^ (29) 

tan (360° 2 r / P + P(=,irr-i(Rj / R 0 ) ) - 90°) 


Now solving for Zr in terms of Y<j> : 


Z r - 


tan 1 j 


360 1 


r; 

y| + R o c DE(?i n" 1 ( R J / R e ) 


t! 


(30) 


t- 90° 


Ps in" 1 (R; / R o ) 



Now eqn. (24) and (30) can be solved simultaneously for Y*. This value 
of is given the name Si , which is the distance from the sleeve impact 

point along Ybb- 

The particle landing position at Si along Ybb can then be 
referenced to the center of the auger by matrix transformation. 

Finally, Rj can be varied from Rtube to Ro to plot a landing path. 
Throughput Prediction 

From landing position after / cycle, an estimate can be made for 
throughput rate, based upon single particle analysis. Using equations 
(24) and (30), a computer program was written (see Appendix) to 
calculate si in the Ybb direction, starting from wall impact. A 12 fence per 
turn auger was overlaid on a plot landing of position (Fig. 16) from 
computer generated data points, which account for the rise of the floor 
(Fig 17). The computer program was only used to generate data from 
Ri = 1 .75 " to 2.50 " since multiple wall reflections were involved, and the 
math model only accounted for one. For 10 particles from figure 16, the 
following data was taken : 


# of particles # of fences jumped % 

2 2 20 

4 3 40 

4 4 40 


Figure 15 shows the distribution comparison with experimental results. 
The average number of fences jumped per stroke may then be 
determined by : 


(.2) (2) + (.4) (3) + (.4) (4) = 3.2 fences / stoke 



Assuming the distribution of sand in the auger is approximately 1 pound 
per turn, and a fence spacing of 12 fences per turn, the distribution of 
sand may be expressed as : 

1 / 12 = .08333333 pounds / fence 

The predicted mass throughput may then be expressed as : 

(3.2 fences / stroke) (.08333333 pounds / fence) = .26667 # / stroke 

This, then represents the mass throughput predicted by the math model. 
At the ideal conditions assumed for the derivation, the mass throughput 
is .26667 pounds / stroke. 

Lunar Application 

The application of lunar constraints to the performance 
characteristics of the auger was accomplished by varying the 
gravitational constant in the path profile program (see Appendix). The 
initial, Earth based model called for a separation velocity of 20.5 ft / sec in 
order to achieve a particle path that was virtually linear. 

When the lunar gravity constant was substituted into the 
program, the required separation velocity was found to the around 4 ft / 
sec. Atmospheric drag forces were considered negligible for the Earth 
model and none existent for the lunar model. 

PROTOTYPE and TEST RESULTS 

The motivation for building a prototype is to determine if the 
concept on paper is feasible on a practical, real world level. If some 
degree of practicality can be demonstrated, then tests can be run and the 
results compared with the predictions of the math model. 

The setup used for prototype testing of the Vertically Reciprocating 
Auger ( VRA ) is described by Figure 13. Specifically, the setup consists 
of the following elements : 



1 ) Auger--The auger supplied to this group is a thro© foot section 
of a steel auger used on earth to drill post holes. It is five 
inches in diameter, with a pitch of five inches (Fig 14). To 
prepare the auger for duty as a VRA, it was sandblasted and 
fitted with a total of twelve "fences" per revolution. These 
fences are simply constructed of thin wooden strips, one inch 
high, attached to the auger service with glue and extending 
from the auger's centerline to its outer edge. Finally, the auger 
is wrapped in a transparent plexiglass sleeve to contain the 
particles as they ascend while allowing observation of particle 
movement. 

2) Test Frame--The purpose of this component is to support the 
VRA and its driving mechanism during testing. A length of 
pipe fastened vertically above the auger runs through its 
hollow center, serving as a guide rod to constrain the augers 
motion to a purely vertical reciprocation. 

3) Driving Mechanism--An air actuator is used to provide the 
force needed to drive the auger vertically. An air actuator was 
chosen because it is much cheaper than an electical or 
hydraulic actuator, and simpler than a purely mechanical drive 
of the "Geneva Mechanism" variety. Since mass transfer 
using the sawtooth wave form model is not expected to be 
frequency dependent, control of the air actuator is 
accomplished using a simple hand valve. The actuator has a 
bore of approximately four inches, and the force produced at 
any given line pressure can be calculated for use in the math 
model. 

Testing consisted of two separate stages, the first in which the 
single particle math model was tested using ordinary steel "BB's", and 
the second in which the mass transport of sand through the VRA was 
observed. Data from these tests is given in Tables 1--4 on the following 


pages. 



The single particle BB test consisted of two parts, the first in which 
ten BB's were lined up along one arbitrary fence within the auger. The 
auger was then displaced once and allowed to sink back in a manner 
consistent with the sawtooth model. The distribution of the ten BB s 
within the auger was then recorded as the number of fences jumped by 
each particle. Approximately 40% of the BB’s ended up at or below their 
starting fence; they jumped backwards and made either negative or zero 
progress up the auger. However, the 60% that made positive progress 
up the auger made a greater relative progress. In other words, the 
positive progress BB's jumped more positive fences than the negative 
progress BB's jumped negative fences. This provides a net positive flow 
of BB particles up the auger. From the overall data, it is determined that 
the average number of fences jumped by the BB particles per stroke of 
the VRA is 0.71 , confirming the theory behind the VRA--that particles will 
indeed ascend the auger on a net basis. 

The second part of the BB particle test consisted of placing a 
single BB at the bottom fence, then counting the number of strokes 
necessary to pass this BB through a total of 24 fences, or two revolutions 
up the auger. If the value of 0.71 fences jumped per stroke, found in the 
first part of the BB test, is correct, then the total number of strokes 
required to move the BB up 24 fences is 24 / 0.71 = 33.8 or 34 strokes. 
The data taken on this second part agrees with this prediction, showing 
an average of 31.2 or 31 strokes required to move the BB 24 fences. 
From the close agreement of the data taken in both parts of the single 
particle model BB test, it is concluded that the tests were uniformly 
performed, with little error introduced by the operators. 

The second stage of testing involved a determination of the 
approximate mass flow rate of sand through the auger. Ordinary 
construction sand, used in the manufacture of concrete, was utilized. It 
was dried thoroughly and purged of any particles larger than 3/4 inch in 
diameter. The bottom of the VRA was filled with this sand, and kept full 
while a number of strokes were applied to distribute sand throughout the 
auger. The depth of sand was brought to a level that was at least 
covering over the top of all fences. A series of tests were then run in 
which a total of twenty strokes were applied while the sand exiting the top 
of the VRA was collected. The data from this test shows that an average 
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of 1 .5 pounds of sand passes through the VRA in twenty strokes, which is 
equivalent to a mass transfer rate of 0.075 pounds per stroke. 

The tendency for the particles to jump backwards, which was 
observed in the BB test, was not nearly so prevalent in the sand test. 
Since the sand particles are much smaller than the BB’s, it is nearly 
impossible to test a small group of them and get an exact distribution as 
was the case with the BB’s. However, observation of the sand as it 
ascends the VRA indicates that backwards jumping of fences is much 
less of a factor in the overall mass transport of the sand than it is in the 
transport of the BB's. 

This difference in behavior can be attributed to the damping 
characteristics of the two different transport mediums. The BB particles 
are quite elastic in their motion, bouncing wildly around inside the VRA 
on pathlengths spanning up to seven fences. When a BB strikes the floor 
of the auger close behind a fence, it has a tendency to bounce up against 
and back off of the rear of that fence, sending it in a backwards direction 
which causes it to become part of the 40% which end up below their 
starting fence. The collisions that the BB makes with the auger surface 
are metal— metal collisions, which incur very little damping. In contrast, 
the sand particles exhibit a great deal of damping. As they proceed 
along their respective pathlengths, each particle comes into contact with 
hundreds of other sand particles. These midair collisions damp out some 
of the energy of the particles, causing them to settle more quickly after 
their initial collision with the ceiling above. In addition, the sand particles 
tend to "coat" the surfaces within the VRA , so that collisions with the floor 
or sleeve are likely to be damped by collisions with other particles 
already on those surfaces. Thus the energy of the system, which was 
initially provided by the thrust of the VRA stroke, is quickly lost through 
damping, and the sand particles settle down without bouncing backward 
down the auger. 

This explanation is affirmed by observation of a test with many BB 
particles. If the bottom of the VRA is filled with, say, 500 BB's, then the 
motion of the BB's through the VRA becomes much more like that of the 
sand, that is to say, the BB particles tend to move as a group up the 
auger with much less of a tendency to bounce backwards. Only when a 
"straggler" BB gets caught below the group, without damping from its 



neighbors, does is tend to exhibit a great deal of negative movement 
throught the VRA. When this happens, its motion becomes like that of the 
BB particles in the first BB experiment in which the distribution of ten BB's 
for one stroke was studied. 

It should be noted that the line pressure used for the BB test was 
35 psi, while the pressure used for the sand test, was 50 psi. An increase 
in pressure was needed to accomodate the weight of the auger full of 
sand. 


CONCLUSIONS and RECOMMENDATIONS 

The VRA device as tested held approximately one pound of sand 
per turn. At the distribution of particles given by the math model, an 
average jump of 3.2 fences/stroke, the mass throughput was predicted to 
be .2667 pounds/stroke. This is a factor of 3.56 times greater than the 
actual throughput, which was found to be .075 pounds/stroke. For the 
test set of twenty strokes, the math model predicts a total throughput of 
5.333 pounds, compared to the actual throughput of 1.5 pounds. 

The reasons for this discrepency are numerous. For example, the 
approximate mass of sand per turn, which was assumed to be one 
pound, is in fact constantly changing and varies greatly. Additionally, the 
math model does not take into account the fact that there is a certain 
amount of zero and negative progress of sand ( though there is certainly 
less with sand than there is with the BB particles ). Damping effects, 
which lower the energy of the particles and cause them to settle quickly, 
are not considered in the math model. In short, the math model 
represents the idealized case; it is not surprising that "real world" 
performance of the VRA is somewhat below predicted levels. Still, the 
basic theory behind the VRA-that there will be a net positive flow of 
particles up the auger--is confirmed. 

The goal of this investigation was twofold: first, to prove the theory 
behind the VRA; second, to propose a mathematical model which can 
predict throughput to some reasonable degree of accuracy. Both goals 
have been met, with the mathematical model differing from the real world 
case by a factor of only 3.56. Based on our experimental results, then, a 



correction factor of 1/3.56 = .281 for the math model prediction is 
suggested. 

The actuator method of VRA drive, used in conjunction with the 
modified sawtooth math model presented in this report, turned out to be a 
simple and effective combination, in terms of both analysis and 
implementation. The setup has, apart from the control valves, only one 
moving part. It is not frequency dependent, and thus requires neither 
tuning nor precise control. Finally, it is compact and, due to its paucity of 
moving parts, quite efficient ( nearly as efficient as the air compressor 
used to supply the air actuator ). A system similar to this one, then, would 
be a good choice for a lunar application. 

The following recommendations are made concerning the math 
model and auger design : 

- Modify the given math model to predict throughput rate for a 
sine wave input, and predict an optimal resonant frequency. 

-.Predict throughput for and test a square wave input that allows 
two particle jumps each cycle. This will require auger throws 
to be greater than the pitch, depending upon auger design. 

- Consider the use of reflectors above the fences that effectively 
increase the pitch of the ceiling. 

- Consider the use of a multiple helix auger that allows both 
steep pitch and low ceiling heights. 

- Explore the possibilities of optimizing such parameters as 
fence height, actuator throw, actuator pressure, sand particle 
size, feed rate and fence spacing. 
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TABLE 2 

PERCENTAGE OF BB PARTICLES JUMPING EACH PARTICULAR 
FENCE (PERCENTAGE DISTRIBUTION OF PARTICLES) 


PERCENTAGE OE PARTICL.ES 
JUMPING EACH PARTICULAR 
FENCE 


TOTAL # PARTICLES JUMPING 
EACH FENCE FOR 10 STROKES 


TOTAL OF 10 BB PARTICLES 


l 00 


— 


-3 

-2 

-i 

1 

FENC 

o 

:es : 

l 

rUMPE 

2 

:d 

3 

A 

L 

5 

l 1 

1 

7 ! 

i 1 

PERCENTAGE 

OF 

PARTICLES 

£ 

i 

1 

1 

8 

1 A 

1 5 

3£ 

1 £ 

1 £ 

1 

2 

1 

1 

0 

1 

1 1 
i 

i 


JUMPING 

EACH 


PARTICULAR 
FENCE 
( V. ) 
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I 


PERCENTAGE GF PARTICLES 
JUMPING EACH PARTICULAR 
FENCE 


X 


FENCES JUMPED 


FENCES 


NUMBER OF 


STROKES 


STROKE 

L 


( .OEM -4 ) + < .01 ) ( -3) 

( . 15 ) ( 0 ) + ( . 35 ) < 1 ) 

( . 02 ) ( 4 ) + < .01 ) ( 5 ) 


( . 08 ) <- 2 ) 
( . 12 ) ( 2 ) 
( 0 )( 6 ) 


+ < . 14 ) <-n 
+ (. 12 ) 0 ) 

+ ( . 01 ) ( 7 ) 


+ 


10 


J 


= 0.71 


i 

FENCES i 


STROKE 


THUS THE AVERAGE NUMBER OF FENCES JUMPED BY THE BB 
PARTICLES PER STROKE OF THE VRA IS EXPERIMENTALLY 
FOUND TO BE 0.71. THIS NUMBER IS LESS THAN ONE, 
INDICATING THAT TO GAIN A NET FLOW OF BB PARTICLES 
UP THE RAMP, AT LEAST TWO STROKES MUST BE APPLIED. 

THIS NUMBER IS LOW DUE TO THE FACT THAT APPROXIMATELY 
40*/. OF THE BB PARTICLES JUMP BACKWARDS, HAVING A 
NEGATIVE EFFECT ON THE OVERAL L NET MASS FLOW OF B B ' S . 
THIS BACKWARDS JUMPING IS MUCH LESS IMPORTANT IN THE 
MASS FLOW OF SAND, DUE TO DAMPING EFFECTS WHICH QUICKLY 
LOWER THE ENERGY LEVEL OF THE SAND PARTICLES AND PREVENT 
EXCESSIVE BOUNCING. 
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TABLE 3 

DATA GATHERED FOR NUMBER OE STROKES NECESSARY TO 
MOVE A SINGLE BB PARTICLE UP A TOTAL OF 24 FENCES 


TRIAL # 

# OF STROKES NEEDED TO JUMP A SINGLE BB 
PARTICLE A TOTAL OF 24 FENCES 

1 

29 | 

j 

2 

27 | 

3 

38 

4 

32 

5 



30 

— j 


THE AVERAGE NUMBER OF STROKES NEEDED TO MOVE A SINGLE BB 
PARTICLE UP A TOTAL OF 24 FENCES IS GIVEN BY: 

29 + 27 + 38 + 32 + 30 

= 31.2 STROKES 

5 


THIS COMPARES NICELY WITH THE NUMBER OF STROKES WHICH WOULD 
BE PREDICTED BY THE EXPERIMENTAL VALUE OF 0.71 FENCES PER 
STROKE, FOUND FROM TABLE 2: 


24 FENCES 


0.71 FENCES PER STROKE 


33.8 STROKES 


FROM THE CLOSE CORRELATION OF THE DATA, IT IS DETERMINED THAT 
THERE IS LITTLE ERROR INTRODUCED BY THE E XPF.R 1 MENT ORS IN THE 
FORM OF NONUNIFORM CONTROL OF THE TEST APPARATUS, OR NONUNIFORh 
MEASUREMENT OF RESULTS. 
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TABLE 4 

DATA FOR THROUGHPUT OF SAND, GATHFRED FOR SEVERAL 
SETS OF 20 STROKES EACH APPLIED TO THE VRA AFTER 
FIRST DISTRIBUTING SAND UNIFORMLY THROUGHOUT AUGER 


TRIAL # 

# STROKES 

! 

WEIGHT OF SAND THROUGHPUT (POUNDS) j 

1 

30 

1.55 1 

| 

2 

20 

1.45 1 

j 

3 

20 

1 .60 

4 

i 

20 

1 .85 

1 

5 

20 

1 .05 


AVERAGE SAND 
THROUGHPUT 


WEIGHT OF SAND THROUGHPUT 
AT EACH TRIAL 


NUMBER OF TRIALS 


1.55 + 1.45 + 1 . 4> + 1.B5 + 1.05 | 

J 


5 


1 


= 1 .5 POUNDS PER 20 STROKES 


= .075 POUNDS PER STROKE 


THE MASS THROUGHPUT OF SAND IS APPROXIMATELY 1.5 POUNDS FOR 
TWENTY STROKES, WHICH COMES TO .075 POUNDS PER STROKE OF 
SAND TRANSPORTED THROUGH THE VRA . 



LANDING POSITION DATA 
FROM 3-D COMPUTER MODE 


R i 

S 1 

1 . 75 

3 . 4 

1 . 8 

3 . 3 

1 . 85 

3 . 2 

1 . 9 

3.11 

1 . 95 

3 . 04 

2 

3 . 0 

2 . 05 

2 . 9 

3 . 1 

2 . 84 

2.15 

2 . 80 

2 . 2 

2 . 77 

2 . 25 

2 . 75 

2 . 3 

2 . 74 

2 . 35 

2.73 

2 . 4 

2 . 75 

2 . 45 

2 . 78 

2 . 5 

3 . 0 



APPENDIX 



PROGRAM SPACE ( I NPUT, OUTPUT) 3~D P/\RT|CL£ P/\T H 

REAL N 1 (500) 

REAL N2 (500) 

REAL X(500) 

REAL Y (500) 

REAL D I SX (500) 

REAL X2.T1 .T2.T3.V1 
INTEGER l,W 

PARAMETER (G=32 . 2 , P= . 4 1 66667 , M= 1 .0 , B=35 . 32* .0 1 71 * 533 ) 

PARAMETER (D=-54 . 68 *. 0 1 71*533, F=375 .0 ,THETA=1 7 . 66 *. 0 1 74533) 

C THIS PART OF THE PROGRAM CALCULATES THE RISE 

C OF THE PARTICLE'S PATH 

X2= (F/M) -G 
T1=((2*P) /X2)**.5 


V1=X2*T1 




PRINT*, 

'VI 

IS 

'.VI 

PRINT*, 

' T1 

IS 

1 *T1 

PRINT*, 

1 X2 

IS 

' ,X2 


C THIS PART OF THE PROGRAM CALCULATES THE FALL IN 

C PARTICLE 

1=0 

DO 10 T2= 0 . 0 , 1 .0 , .0005 
1 = 1 + 1 

N 1 (I) = ( (G*COS (B) *T2**2) /2)+Vl*T2 
N 2 ( I ) = (-G*SIN(B)*T2**2)/2 
IF (N1 (I) .GE. P) GO TO 20 
10 CONTINUE 

20 W=0 

DO 30 W= 1,1,1 

X(W) =N 1 (W) *COS (D) + N2 (W) *C0S (B) 

Y (W) =N1 (W)*SIN (D) + N2(W)*SIN(B) 

PRINT*, 1 X= 1 , X (W) * 1 2 
PRINT*, 1 Y= 1 , Y (W) * 1 2 
D I SX (W) = X(W) /COS (THETA) 

PRINT*, 1 D I SX= ' , D I SX (W) * 1 2 
CONTINUE 
PRINT*, X (1-1) 

PR I NT*. Y (1-1) 

STOP 
END 


30 



OPTION BASE 1 
DIM A(2 , 50) 
REAL Rl,Ro, P 
REAL Asp, Beta 
INTEGER I 
1=0 


3~D Ianoii y<7 
am 


Ro=2 . 5 
P=5 . 0 

FOR Rl=. 9 TO 2.5 STEP .05 


PoSiT/oa/ 


1 = 1+1 

Asp=ASN (Rl/Ro) 

Beta=2 *ATN ( P/ ( PI *R1 ) ) 

FOR Z=. 02 TO 5 STEP .01 

Yr=Ro*COS (Asp) - (Rl/ (TAN (2*PI*Z/P-PI/2+2*Asp) ) ) 

Yp=(Z-P+Rl/ (TAN (Asp) * (TAN (Beta) ) ) ) * (TAN ( 2*ATN (P/ (2*PI*R1) )))*(-!) 

IF ( Yp-Yr) > . 05 THEN GOTO 200 

A (2 , I) =Yp 

A ( 1 , I ) =R1 

PRINT Rl , Yp 

GOTO 210 

NEXT Z 

NEXT Rl 

END 



Progress Report - Week 1 


Linear Auger Group 

This week we familiarized ourselves with the project itself, which is a 
reciprocating auger used in drilling and removal of debris at lunar excavating sites. 

1. Etheridge, Mark - Reviewed previous reports. 

2. Fair, Robert - Organized meetings and became familiar with the Apollo system. 

3. Morgan, Scott - Became familiar with VersaCad system and looked through 
Lunar Bases book for references. 

4. Pearson, Brent - Contacted suppliers for a plastic sleeve. 

5. Weldi, Kevin - Reviewed previous reports and became familiar with the Apollo 
system. 

6. Woodrough, Steve - Reviewed previous reports. 

Due to inclement weather, the group met Tuesday at 3 pm to discuss the project 
and scheduling. 



MEMORANDUM 


TO : Dr. Brazell 

FROM : Linear Auger Group (C) 
DATE : 1-21-87 

SUBJECT : Week #3 Progress Report 


On Thursday 1-14-88 at 8:00 PM the group met with Dr. Brazell and 
discussed the different parameters, constraints and objective concerning the 
auger design. 

On Wednesday 1-20-88 at 12:00 PM the group met again to begin work 
on the mathematical models and discuss variation of fences and possible 
reflectors. Discussed methods of welding and attaching fences/reflectors. 

Individually : 

Robert Fain: Worked on Versacad. 

Scott Morgan: Worked on Versacad. 

Kevin Weldi: Researched Project more 
Steve Woodrough: Worked on dynamics and math model 
Brent Pearson: Learned word processor use (PC write) and how to get laser 
print. 

Mark Etheridge: Worked on math model 


MEMORANDUM 


TO : Dr. Brazell 

FROM : Linear Auger Group (C) 
DATE : 1-28-87 

SUBJECT : Week #4 Progress Report 


On Thursday 1-21-88 at 8:30 the group met in the Design Lab for 2.5 hours 
for discussion of problem statement, math model, actuating device and 
division of responsibilities. 

Arranged to have auger surfaces polished and edges ground. 

Group met 1-27-88 at 12:00 PM to discuss integration of system elements. 
The group was then divided into four sections. 

Individually : 

Robert Fain: Worked on Versacad. 

Scott Morgan: Math model and frequency determination. 

Kevin Weldi: Power Requirements. 

Steve Woodrough: Spoke with Vendors, GTRI, and met with Brent. 

Brent Pearson: Worked on drive mechanism. 

Mark Etheridge: Math model and frequency determination. 



MEMORANDUM 


TO : Dr. Brazell 

FROM : Linear Auger Group (C) 
DATE : 2-04-87 

SUBJECT : Week #5 Progress Report 


Group C met in the lounge of the Coon Building at 12:00, Wednesday Feb. 13. 
The math model, as well as the relationship between math model and test 
setup, were considered. 

The decision has been made to design the test setup using an air actuator, if 
possible. This device would provide an input more closely approximated by 
a square rather than the sine wave which had originally been anticipated. The 
math model therefore, has been reviewed and found in need of 
modification. A new math model is being developed at this time. 


Individually : 

Robert Fain: Worked on Versacad. 

Scott Morgan: Math model and frequency determination. 

Kevin Weldi: Worked with Data-Grapher and Graphics Editor. 

Steve Woodrough: Spoke with Vendors, GTRI, and met with Brent. 

Brent Pearson: Drew plans for test-setup, determined materials needed. 
Mark Etheridge: Math model and frequency determination. 



MEMORANDUM 


TO : Dr. Brazell 

FROM : Linear Auger Group (C) 
DATE : 2-11-87 

SUBJECT : Week #6 Progress Report 


Progress this week centered around preparation for the mid-term presentation. 
Scott Morgan was chosen to give the presentation. The various sketches, 
drawings, and graphs were compiled and put on transparencies for class 
viewing. Of the several possible math models, one was chosen for use in 
testing procedures, mainly due to the simplicity of its implementation. Vendors 
have been contacted and plans made for fabrication to begin this weekend. 

Individually : 

Robert Fain: Worked on Computer programs and Versacad drawings. 

Scott Morgan: Prepared for and gave presentation. 

Kevin Weldi: Worked on Datagrapher and Versacad drawing. 

Steve Woodrough: Helped with Math model and contacted vendors. 

Brent Pearson: Worked on Versacad and communicated with vendors. 

Mark Etheridge: Continued work on Math model and Computer program. 



MEMORANDUM 


TO : Dr. Brazell 

FROM : Linear Auger Group (C) 
DATE : 2-18-87 

SUBJECT : Week #7 Progress Report 


This weeks focus was on the construction of the test stand and methods 
of coupling the auger and actuator to allow two degrees of freedom to 
prevent binding and reduce vibration. 

The group also discussed the final report setup and reviewed the patent 

forms. 

The drawings used in the oral presentation are being translated to 
computer for further analysis. 

Considering computer controlled solenoid valve for testing of a the 
math model. 

Individually : 

Robert Fain: Working on Computer controlled actuator design. 

Scott Morgan: Worked on Math model. 

Kevin Weldi: Revising drawings. 

Steve Woodrough: Beginning work on written report. 

Brent Pearson: Built test stand. 

Mark Etheridge: Continued work on Math model and Computer program. 



MEMORANDUM 


TO : Dr. Brazell 

FROM : Linear Auger Group (C) 
DATE : 2-25-88 

SUBJECT : Week #8 Progress Report 


This week we found the fabrication of the test setup slowed due to 
problems with the air actuator. Line losses have caused the speed with which 
the actuator throws to be too slow. A less restrictive air delivery system, of an 
actuator requiring a smaller volume of air, must be investigated. The fences 
for the auger have been purchased and will be applied immediately. Final 
fabrication is planned to be completed this weekend, so that testing of the 

model can begin. . 

Work continued on the math model, with problems arising in the 
manipulation of 3-D vectors. An Alternative to the purely mathematical vector 
method of tracing the path of a particle as it proceeds up the auger, suggested 
by Mr. Brazell, is to use wire of thread to physically trace the path on the auger. 

Individually : 

Steve Woodrough: Continued design on the test setup and contacted venders. 
Scott Morgan: Continued work on the math model and physical representation 
using stiff wire. 

Mark Etheridge: Worked on math model with Scott. 

Robert Fain: Continued work on drawings, computer program and test setup. 
Kevin Weldi: Overhauled and tested actuator and worked on design setup. 
Brent Pearson: Continued design and fabrication of test setup, purchased 
materials. 
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DISCLOSURE OF INVENTION 
SUPPORTING INFORMATION 


Setting aside your personal interest, what do you see as the greatest 
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Alternate Technology and Competition: 
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the purpose of the invention. 
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use of these alternate technologies. 
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